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We evaluate the energy momentum tensor of a massive quark as it moves through 
^ I an A/ = 4 SYM quark gluon plasma at constant velocity. We find that in the 

Q"^ ■ near-quark region, where the dynamics is expected to be dominated by dissipative 

CN ■ behavior, the energy density may be quantitatively characterized by a transient at 

> ; 

' velocities above the speed of sound of the plasma. 
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5 ■ I- INTRODUCTION AND SUMMARY 

Q-i! The Relativistic Heavy Ion Collider (RHIC) at the Brookhaven national laboratory 

■ 3,4] (see for a review) probes the dynamics of an expanding quark gluon plasma (QGP) 

> '• 

•'-j . at temperatures shghtly above the QCD deconfinement temperature Tqcd ~ 170 MeV. One 
r> ' 

^ I of the prominent features of the QGP at RHIC is the quenching (lack of energy) of hadronic 
jets emitted from the plasma. It is understood that the plasma acts as a dissipative medium 
for the partons moving through it. The exact mechanism for the energy loss of the partons 
is unclear partly because the coupling constant gvM is quite large and certain features of 
the QGP are not accessible perturbatively, and partly due to the difficulty in modeling the 
interaction of a parton with the surrounding quark gluon plasma. 

One approach which overcomes these difficulties is to consider, instead, a quark in a 
strongly coupled A/" = 4 large SU (N) SYM plasma. While there are several differences 
between the former and a QCD plasma, one of the lessons that has been learned over the past 
year is that many features of an A/" = 4 SYM plasma roughly coincide with those of the QGP 
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at RHIC. A comparison includes the energy density, the viscosity to entropy density ratio 
Q], the friction coefficient for the moving quark and the jet quenching parameter flol.[ll| 
One may extend the analysis to include theories which break conformal invariance or A/" = 4 
supersymmetry. For example 12| and 13| study the viscosity to entropy density ratio in a 
non conformal theory and a theory with massive quarks in the fundamental representation. 
The diffusion coefficient for various deformations of the superconformal Af = 4 theory have 



been calculated and analyzed in 



Q, H, Q, 0, Q, 
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Similar investigations of 
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301 ] ■ Some 



the jet quenching parameter have been carried out in 
other features of the quark gluo n plasma in AA = 4 SYM and related theories have been 



studied in 
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This is not to say that M = 4 SYM theory is a good approximation to QCD. While there 
is qualitative agreement between the two, on the quantitative level there is the expected 



mismatch (see 



11| and [42| for some discussions). Here, we shall adopt a somewhat con- 



servative approach and try to understand the characteristics of the energy deposition of a 
moving quark in an A/" = 4 SYM plasma as a ffist step towards understanding quark-QGP 
dynamics. While being motivated by the comparisons between A/" = 4 SYM and the results 
of the RHIC experiment listed above, whether our analysis is actually related to QCD or 
not remains somewhat speculative. 

To study the energy associated with a quark we shall evaluate the near-quark energy 
momentum tensor as the quark travels through a strongly coupled QGP at constant velocity. 
This may be done by appeahng to the AdS/CFT duality [43|. There, one may evaluate 
the energy momentum tensor of a heavy quark in A/" = 4 SYM at finite temperature by 
considering metric fluctuations due to a string han ging down from the boundary of an AdS- 



39, 



44 



45| 



Schwarzschild (AdS-SS) background geometry 

We are interested in the temperature dependent contribution of the quark motion to the 
energy momentum tensor. Apart from this subleading term, the energy momentum tensor 
will also be influenced by the leading zero temperature contribution of the moving quark, ^ 
and by the QGP itself. Our main result is the real space value of the energy momentum 



^ The value of the leading zero temperature contribution of the moving quark to the energy momentum 
tensor may be obtained up to an overall constant by imposing conformal and Poincare symmetry 39|. 
We reproduce this result in equation ([5T|) . 
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tensor {T^^) (the 'd' is a reminder that the energy momentum tensor of the conformal 

d 

plasma and the leading contribution from the quark motion at zero temperature have been 
subtracted) which is valid in the vicinity of the quark 
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have been omitted and contributions with 



In equations ([T]) C(T^) corrections to {T^u) 

d 

delta-function support where ignored. 

Our notation is such that x_ = x — f t is the comoving coordinate of the quark, T is the 
temperature of the quark gluon plasma, v the velocity of the quark (in the x direction), r 
is the direction transverse to the quark motion, 9 is the polar angle and A is the large N 't 
Hooft coupling constant. 

An analysis of (T^jy) was initiated in 39| where the large distance behavior of the energy 
momentum tensor was discussed. There it was shown that if the quark moves faster than 
the speed of sound, = ^ Mach cone will be created far from the moving quark. 



or the 



461] when 



This fits nicely with a hydrodynamic model [461] (valid at large distances as wel ^ 
sound waves carried by the plasma as a result of the quark motion. As argued in 
one gets closer to the moving quark he or she will reach a regime where hydrodynamics 
becomes non-linear and turbulent effects are expected. In the immediate vicinity of the 
quark there is a region where the hydrodynamic approximation breaks down completely and 
strong dissipative effects dominate the dynamics. This is precisely the region described by 
equations ([1]). It is nice to see that the AdS/CFT duality provides for a method of probing 
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the near-quark surroundings which is not accessible via standard constructions due to strong 
dissipative behavior (see 46|]) and the lack of a model for the relevant interactions. 

Focusing on the energy density {Ttt{x-,r)) we find that it exhibits some velocity de- 

d 

pendent directional features (see figure [1]) which do not appear for the other components of 
the energy momentum tensor. As long as the velocity is relatively low, < < 5/13, the 
energy density has a vanishing gradient in the x_ direction along 

tan uji = [ — = ^ 1 — V ). 

\x-J -20 + 44t;2 ^ ^ 

In this range of velocities no exceptional behavior is observed even when the quark passes 
the speed of sound = 1/3. As mentioned earlier, this should be contrasted with the Mach 
cone which is formed far from the quark as it passes the speed of sound of the conformal 
plasma 



39 



46|. 



At velocities squared above 5/13, the subleading dissipative contribution to the energy 
density begins developing lobe-like features specified by a non trivial, vanishing gradient in 
the r direction at an angle of 

^ 5 - 13v^ + 8v^ 



tan UJ2 = — 



x_ J —5 + ISv"^ 

As the velocity squared of the quark is increased to 5/11 one finds that there are two 



contributions to {Tu) 
a new low ene 
at an angle of 



the previous contribution now deformed to a lobe of high energy, and 

d 

a new low energy contribution right behind the quark. ^ This new contribution is directed 



2 I r \ 5 + + VTEVTE^^bOv^T^lv^ , 2^ 
tan uJs = { — = ( 1 — f 



As the velocity is increased further, the energy deficit behind the quark dominates and the 
contribution of the lobe becomes smaller until it vanishes at a velocity squared of 5/8. We 
encourage the reader to have a look at |i48i] where this effect can be observed as a continuous 
function of the velocity. 

This behavior may suggest that there are two mechanisms that come into play. One at 
relatively low velocities v"^ < 5/13, responsible for a region of high energy, and one at high 



^ The terms high and low energy are, of course, relative to the leading zero temperature contribution of 
the quark and the contribution of the plasma which have been subtracted. We would like to thank the 
authors of 47| for clearing up a sign mixup in this context. 




FIG. 1: Contour plots of the leading dissipative contribution to the energy density 
J- {Ttt{x-, r)) due to the motion of a quark moving at constant velocity v. The horizon- 
tal axis corresponds to the comoving coordinate of the quark, X- = x — vt, while the vertical one 
corresponds to the direction transverse to the quark motion, r. Large positive values are shaded 
in red and marked with an 'H' while lower values correspond to blue shading and are marked 'L'. 



See 



4811 for an animation. 



velocities v"^ > 5/8, responsible for a region of energy depletion, with some overlap. It is 
interesting to note that the overlap region is above the Mach velocity for the quark, and that 
there is no remnant of the Mach cone discussed in 39|].^ A different possibility is that in 



^ Note however, that when the Mach cone is generated there is a buildup of energy in its vicinity (see 
figure 2, and equation (80) of [39|) which may extend all the way down to the near quark region at the 
intermediate velocity range 5/11 < < 5/8. 
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figured] we are seeing an interference pattern of some sort. Though then it is not clear what 
sort of waves are propagating in the plasma since sound waves are not expected at such high 
momentum scales and the gluon field strength does not exhibit lobe-like behavior at short 



distances 4l| . It is perhaps worth mentioning that a region of energy deficiency behind the 



quark is characteristic of an LPM (Landau- Pomeranchuk-Migdal) effect 49|, l50|] where the 
parton motion in the plasma is modeled by inelastic collisions with fixed scatterers. The 
literature contains many more interesting discussions regarding the mechanism for energy 
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51 



52 



53|) though we could not find quantitative 



loss in a QGP (see for example 
correlations between the descriptions there and the behavior described above. 

While contemplating over these issues, the authors of [47|] brought to my attention that 
their work has some overlap with this paper. I thank them for sharing their results with me 
prior to publication. 

II. EVALUATION OF {T^^). 

We wish to evaluate the expectation value of the energy momentum tensor in response to 
the motion of a quark in a thermal plasma in A/" = 4 SYM. Using the AdS/CFT dictionary 
this is given by the response of the AdS-SS metric to the motion of a string with one endpoint 
located on the AdS boundary, tracing the worldline of the quark. 

Ignoring for the moment the effect of the string on the spacetime geometry, one can find 
the string's profile by minimizing the Nambu Goto action 

1 



Sng = J V-wdrda 

in an AdS-SS background. Hence, w is the induced metric Wa/s = G^/^udaX'^d/sX" with G^^^ 
the metric of an AdS-SS black hole; G^^^ may be read off of the line element 

ds"^ = — —q(z)dt^ + dx^ + > dx^ H 

with g{z) = 1 — (y'^ ■ The coordinate z runs from at the asymptotically AdS boundary to 
zq at the black hole horizon. We have set the value of the dilaton to its trivial background 
value since it decouples from the equation of motion at the linearized order that we are 
considering. 
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The equations of motion for the embedding functions of the string X^{a, r) are supple- 
mented by boundary conditions which require that the string endpoint moves at constant 
velocity v, and that no energy is fed to the string from the horizon. Using a static gauge 



and an ansatz X^{a, r) = (r, ^(cr) + vt, 0, 0, a) one finds 



BE 



e(^) = ^(^ln^ + 2arctan^j. (2) 

The differential equation describing the combined dynamics of the string and the target 
space metric (within the supergravity approximation) can be found by minimizing the action 

S = Seh + Sng (3) 

where 



(4) 



Since we will be looking for the expectation value of the energy momentum tensor, it will be 
enough to consider the linearized equations of motion for the metric fluctuations h, defined 
through G = G^^^ + h, in the presence of the string whose profile is determined by equation 

(ED. 

The linearized equations for h resulting from minimizing ([3]) may be written in condensed 
form as 

V^^P^hp^ = T^"" (5) 

with T>'^'^P" a second order linear differential operator. The boundary conditions for the 
metric fluctuations are that they vanish at the AdS boundary (implying that there is no 
source term for the metric) and that there are no outgoing modes coming from the horizon. 
Once a solution, /i, to equation ([5]) is obtained the energy momentum tensor in the SYM 



theory (T) = (Tpiasma) + {Tquark) may be determined through [39| 



{T,uark) = ^\jY^Q^ (6) 

where Q^y is obtained from h = . . . + Qz^ + . . . through a small z expansion. {Tpiasma) = 
T^]\j-2rp4: (li(iqQnal(3 1 1, 1) + 0{N'^) is the contribution of the thermal plasma to the energy 



momentum tensor 



5J] . Since we are interested in the dissipative effects of the quark motion, 



we subtract this background value and so, for our purposes it may be ignored. 
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In equation ([6]) one should use the relations 

nL^ 1 



2G5 



and 



a' 



(7) 



Apart from that, one should also make sure that all lower order terms in the expansion of h 
in small values of z correspond to contact terms (and therefore may be renormalized away^). 

In 39] it was shown that by choosing a gauge such that /is^ = and working in cylindrical 
coordinates, the coupled equations of motion for the various components of the metric tensor 
decouple into three sets of second order equations and two sets of first order constraints. 
Explicitly, we Fourier transform the metric fluctuations, 

h{k,z)= I V(a;,z)e-'('^-"-+'=^"2+'=3-^)d2A;, 



and define 



hfMuik--, k±, z) 



1 



Hfj,u{k-, k±, z), 



(8) 



a'Vl - 

where we have also rotated h to a cylindrical coordinate system. By considering the combi 
nations 

k 



A 




and 



E 



E2 



^ (^01 + f^^02) 
I (-f^ll + -f^22 + 



(9) 
(10) 

\ 



Hu 



+ H22 -^ + Q 



k- 



+ 2H» - 12 



•'#i^i2 



^ A holographic renormahzation prescription was presented in 55, 5^, see also [s?! for some recent applica- 
tions. If one is careful when taking the limit where the AdS radial coordinate reaches the asymptotically 
AdS boundary then a subtraction of divergent contact terms will give the same results. See for example 



58, 



59, 



60, 



61|. 
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one finds that (I5l) reduces to 



Id't+{--I+^Kx]d, + Vxiz)]X 
z g 



(12) 
(13) 



with X = E,D or A and g' = dzg{z). I is the identity matrix in appropriate dimensions. 
The matrices Kx, Vx, Kx and Vx and the vectors Sx and Sx can be found in 39|| and 
have been reproduced in appendix |A] for completeness. An important feature of equation 
( fT2l) that we will use is that in the zq ^ 00 limit^, ^Kx 0. 

Due to the linear relations between E, D, A and H, equations ([9]) to ( fTTl) . the fourth order 
contribution to a small z expansion of hf^^, can be inferred from the fourth order coefficient 
in an expansion of the fields X. Defining X = Nxz + Oxz^ — \Pxz^ + Qxz^ + 0{z^) (recall 
that h{^) = 0), we find that Nx = Ox = and that the first order constraints f|T3|) imply 
that 



Pa — Pdi — Pd2 — Pe2 — ~Pez 



Pe, 



1 + v\ P, 



Ei 



(14) 
(15) 



Inverting the relations between the A, D, E variables and the H variables, equations ( 
to ( fTTl) . one may obtain the third order term P^^, in the small z expansion of h, h^y 

-\P,uZ^ + Q^.z^ + 0{z^), 



P 



/4 + 2f2 -6f \ 

-Qv 2 + 4v^ 

2 - 2v^ 

\ 2 - 2v^j 

Since P is independent of the momenta, any divergences that it may cause will be contact 
terms. This allows us to use equation ([6]). 



^ One should use the dimcnsionless quantity fczo instead of zq (with k some scale with units of inverse 
length) when considering an expansion at "large" values of zq- At the moment, we shall refrain from 
writing this out explicitly to avoid cluttering the notation. 
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The constraint equations (fT3!l also relate the values of Qx to one another. One finds 

1 



Qd2 
Qe2 



Qdi - 



Q 



1 
2 



-1 + 3v' 



Qe, + 



1 



jQEi, 

3ik^v{l + v^) 



(16a) 
(16b) 
(16c) 



Therefore, we need only find Qei, Qdi and Qa to obtain all the fourth order terms in an 
expansion of X. Then, as discussed in the previous paragraph, one may invert equations ([9]) 
to flTTl) to obtain the fourth order term Q^,^ in an expansion of hf^^. Plugging it into equation 

we will obtain the Fourier transformed energy momentum tensor we are looking for. To 
make our notation slightly more precise, we should have used Q^u for the fourth order term 
in the expansion of h instead of Q^u which was defined below equation ([6]) to be its real 
space counterpart. However, since in the rest of this section we will be working in Fourier 
space, we choose an unhatted notation which is less cumbersome. 

To proceed, we note that the equations of motion ( fT2l) have the property that in the 
zq —>■ oo limit, the kinetic terms become diagonal and correspond to the kinetic terms 
obtained for a scalar field in an AdS background. This implies that in the 2:0 ^ oo limit it 
is more useful to change to a basis where the potential is diagonal as well. Our approach to 
solve f|T2l) and f|T3l) is to first redefine the variables D and E so that in the Zq 00 limit we 
shall get the AdS equations of motion. To do this, we multiply the matrix equation ( JT2l) for 
X = D from the left by 

„2 



d 



D 



— (a 



V2[a 
2 , 1^ 



1 



+ 1 



[a 



r 



1 



(17) 



and for X = E hy 

/|(a2 + l)(2a2_i) _4^2(2tt2_i) _4 ^ 1) ^2a2 _ 1) 1 (1 _ 2a2)2\ 
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a' + l) 
3 {a" + 1) 
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2 I ^2 



1 

a2+l 

-4a2 



1 

fa' + 1) 



|(l-2a2) 







I (a' 



2^2) 
+ 1) / 



(18) 



In equations (fTTIl and (ITSi) we have defined a = with k'^ = k'j_ + {1 — v'^)k'^_. 

This will diagonalize the equations of motion in the zq ^ 00 01 g ^ 1, g' Q limit. 
Actually the equation of motion for E can not be diagonalized. Instead we have brought it 
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into Jordan form. The precise equations of motion one obtains are rather messy, and we shall 
not write them out explicitly. What is important, is that in our new basis we are assured 
that the "off diagonal" terms are, by construction, at least of order Zq"^. Therefore, we 
may first solve the diagonal equations and then use perturbation theory to find corrections. 
These will be at least of order Zq^ and will be neglected. Changing to the more convenient 
variables Z = kz, Zq = kzQ, we can concentrate on the following three equations 
.2^f 3 4Z3 Y Z^,{Zl-Z^{a^ + m Ze-^-^(^A) 

3 Z-a- \ Zl{Zl-Z^{a- + l)) \ _ v/2Ze-^-^(^/^) 

(21) 



z S{Z^-Z^) J V {z^-z^f 

Ze-'^-^^^l~^) ( (2^2 _ 1) ((q;2 + 1) ^2 _ ^ 2) 2t;2z4 {o? + 1) (2^2 _ i) 



where 



k-'^Da = V2 [Di [a^ + 1) - D2a^) , (22) 
= -^^2a^ + ^i?4 (2a^ - 1) + ^i^i (a^ + l) - ^E, {a^ + l) , (23) 
k-^A' = A. (24) 

The Da and equations are the first diagonal components of f[T^ after multiplying it from 
the left by di), equation fITTI) . or by c/^;, equation (ITSl) . respectively. 

The y4' equation of motion is identical to the equation of motion for the dilaton in an 
AdS-SS black hole background as it is sourced by the trailing string [38|, l39|, |4l|, |62|. In 



4l| the large momentum asymptotics of equation (fT9|) were found and the near quark value 
of the field strength, ~ TrF'^, was evaluated. The WKB method was used to solve the 
homogeneous equation at large momenta; the solution to the non homogeneous equation 
then follows by constructing the Greens function from the homogeneous solutions. 

Alternatively, to obtain the large Zq asymptotics of A' one may expand the equation of 
motion (fT9l) and the function A'{Z) in an inverse power series in Zq. The unique normalizable 
solution to order Zq'^ is given by 

A\Z) = ^ {Z{h{Z) - L,{Z)) + 2L,{Z)) - (25) 



2k' ' ^ ' ' " " " 3Z 
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where I2 is a modified Bessel function of the first kind and Lj are modified Struve functions 
of order i. A caveat in this approach is that the boundary conditions are imposed at Zq 
which is the perturbative parameter we are expanding in. To properly impose the boundary 
conditions on generic solutions one should resum the whole power series. Nonetheless, the 
solution obtained in this somewhat hand-waving manner coincides with the one obtained 
using the WKB method of |41|]. In both cases, one finds that for 

A' = ... + Qa'Z^ + ..., (26) 

Qa' is given by 

kQA' = ^n-^ + 0{Z^'). (27a) 



16 3Z 



In order to find the 0{Z^) coefficients Qx in the small Z expansion of X = Da and 
X = Ea (in analogy with equation fl26|) ). one needs to solve the equations of motion for Da 
and Ea, (120|) and (12T|) . Following the somewhat heuristic approach used to obtain fl27al) . 
one finds that an expansion of Da-, Ea and their equations of motion in an inverse power 
series in Zq leads to equations of motion which are identical to equation (fT9l) for A up to a 
different normalization of the source terms. Using equations (125|) and fl27al) . one may easily 
read off the coefficients 

Qd^ = V2Qa', (27b) 

_ / (2a^-l)(aV-l) + (^^ + 2)) \ 
Qe. - ^ ^ia^) ) (27c) 



where terms of order 0{Zq^) have been neglected. A more careful analysis following 41 | 
which leads to the same result has been left to appendix [Bl There we also find that the 
regime of validity of equations flTTI) is given by large Zq and a < Z^l"^ . 
Putting together ([16]), ([22]), ([23]) and ([27D we find 



Qd, =Qa + (29) 

Akzq 

Qe, = ^ ((2 + v') - (1 - v')a') Qa + ^, (^(1 + 5.;^) - \a\l - v')^ . (30) 

As discussed earlier, we may now find the contribution Q^^ to the metric fluctuations, 
i^i,, due to the string by inverting the relations ([9]), (fTOjl and (fTTl) . Plugging the resulting 
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Q^y into ([6]) one finds that the contribution of the quark motion to the Fourier transform of 
the energy momentum tensor J {fguark + ^ {{^)) is given by 



k 

24 



/_(2 + t;2_(i_^2 



where 



3v 



3v- 
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-(1 
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P 
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- 2v^) + 

{l-v^)k^k 





{l-»;2)fc2 

p 



(l~v'^)vk±k- 
{l~v^)k^k. 

P 









1 



V 



/ 



gl-L ^ )^ 2k- 



k-k±{l—v'^)v 
9P 



2fc_ 



gv- " fc2 
9p 





2k. 



V 



2\ k_-k±_ 
fc2 



Pk, 



2\ k^k^ 



fc2 



1 1 fc^ 

— (5 - llt;^) + -(1 - v'^W^ 



Pk. 
Pkj 



1 k"^ 
-X5v'^ - 11) + -(1 
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A;2 



1 1 ^2 

2^ ^9^ ' k^ 
18^ ^ 



o\ 




PeJ 
(31) 



(32) 
(33) 
(34) 
(35) 



It is satisfying to see that the leading T = contribution due to the quark motion is exactfy 
what was predicted in 39| based on conformal invariance and Poincare symmetry. Fourier 
transforming the leading dissipative contribution to position space, we obtain equation ([1]). 
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APPENDIX A: THE EQUATIONS OF MOTION FOR THE ENERGY 

MOMENTUM TENSOR 



In order to decouple the components of equation ([5]), we have defined in equations ([9]) 
through ( fTTi) the variables A, D and E. The resulting equations of motion were given in ( fT2|) 
and f|T3|) . The matrices Kx, Vx, Kx and Vx and the vectors Sx and Sx in these equations 
are given by 

2 



K 



A 



K 



D 




Va 



9^ 



9^ 



k 



-9 9v'm 

-1 -Mf)' 



Sa 



S 



D 



(Al) 
(A2) 



/§ o\ 

I 



\0 1/ 



( 



1 + 7 
1 

-l + t;2- ^ 



-1 + 3 I X 



Vf. 



2 / 



6f 






V 2^7 






-12t;2 f ^ 



2^ 

-2^ 






9 
— I 



(A3) 



(A4) 



D 



1 , 4=1 








1 









2 



-S9' 



9 



-3gg' 18v'%g' 3{3v^'$ + g)g' 
\2k'^gz -l2kWz -2 {3v^kl - gk'^) z 2k^gzJ 



(A5) 



(A6) 



(A7) 
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APPENDIX B: WKB APPROXIMATION 

In this section we would like to solve the equations of motion (!20l) and (12T]) with boundary 
conditions such that Da{0) = Ea{0) = and that Da and Ea have only ingoing modes at 
the horizon. In effect, we are not interested in the full solutions to the equations of motion, 
but rather in the fourth order terms in a series expansion of Da and Ea in small Z. 



This problem has been addressed in 41j for the related equation f|T9|) . As discussed in 
the text, the latter equation of motion has been solved by first obtaining the homogeneous 
solutions and then using them to construct the appropriate Greens function. Here we shall 



follow the same method, using some of the arguments made in 41 1 



First, we shall find the homogeneous solutions to equations ( l20i) and ( l2Ti) by bringing 



them to Schrodinger form, ijj" + Vijj = 0. This may be done by defining il^x = \l ^^^^p-X 



with X = EdjDb, and 



bDAZ) = 9iZr\ (Bla) 
bESZ) = giZr'-H'^'^'), (Bib) 



VeAz)-v. 9zMzy ' (^2^) 



and dividing the resulting equation by ^Jj^^^g{Z). The Schrodinger equations obtained 
this way are characterized by the potentials 

VnSz)-v. zMz?^ ^^'^^ 

QAa" {a" + 1)' 



where 



In equation (1B2cl) . Vq is the Schrodinger potential corresponding to the equation of motion 
for a massless scalar field in an AdS-SS background. 

In the notation of equations flB2l) . it is clear that the Vx's differ from Vq only in their 
near horizon asymptotics. In [41] it was shown that as long as the Schrodinger potential 
describing the asymptotically AdS system has a parametrically long flat region for which 



4 / \ 6 
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V{Z) ~ —1 followed by a V{Z) > region near the horizon, then the 0{Z^) coefficient Qx 
in a small Z expansion of X is given by the Z independent terms in 



iD o Jz X 

2/3 

where exponentially suppressed terms, e~^o ^ should be consistently neglected. J{x) is the 
non homogeneous term in the equation of motion for X. In our case (equations (1201) and 
([2ID) we find 

JdSZ) = V2MZ), (B4) 

^^"(^^ = [ 3(^^TT^ 3 {z^ - z^) ) 

where 

Jq{Z) = ie-^'-^^^/'^^giZ) (B6) 
k 

is the source term for the equation of motion for the A' component of the metric (which 
coincides with the equation of motion of the dilaton). b{x) is given in equations (IBll) and 
K2{x) is a modified Bessel function of the second kind. 

Clearly, as long as a < Zq^^ the potentials (]B2al) and (lB2bl) satisfy the flatness and 
positivity requirements discussed above. Therefore, one may use formula (1B3P to obtain the 



fourth order coefficients of Da-, Ea and A'. The A' case has been discussed in [4l|, where 
it was shown that as long as Zq < o? one may expand the integral in (1B3P in a power 
series in Zq owing to the exponential decay of the modified Bessel function K2{x). One 
finds kQA' = jqT^ — + 0{Z^^). To obtain the other two solutions, we note that up to 



order 0{Zq ), the expressions Joay^Da ^'^^ ^Eay^Ea differ from Ja'v^ = Joy 9{Z) by a 
constant multiplicative term. This leads us to equations ([2j 
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